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Abstract 

We compute the next-to-next-to-leading order spin-orbit contributions in the total energy flux 
emitted in gravitational waves by compact binary systems. Such contributions correspond to the 
post-Newtonian order 3.5PN for maximally spinning compact objects. Continuing our recent work 
on the next-to-next-to-leading spin-orbit terms at 3.5PN order in the equations of motion, we 
obtain the spin-orbit terms in the multipole moments of the compact binary system up to the 
same order within the multipolar post-Newtonian wave generation formalism. Our calculation of 
the multipole moments is valid for general orbits and in an arbitrary frame; the moments are then 
reduced to the center-of-mass frame and the resulting energy flux is specialized to quasi-circular 
orbits. The test-mass limit of our final result for the flux agrees with the already known Kerr black 
hole perturbation limit. Furthermore the various multipole moments of the compact binary reduce 
in the one-body case to those of a single boosted Kerr black hole. We briefly discuss the implications 
of our result for the gravitational-wave flux in terms of the binary's phase evolution, and address 
its importance for the future detection and parameter estimation of signals in gravitational wave 
detectors. 
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I. INTRODUCTION 



Previous works [l|, |2|3 have derived the spin-orbit effects in the equations of motion of 
compact binary systems (made of two black holes or neutron stars) at the next-to-next-to- 
leading order beyond the dominant level. Such an order corresponds to the post-Newtonian 
(PN) order 3.5PN ~ in the case of maximally spinning compact objects, i.e. 2PN 
beyond the leading spin-orbit effect at order 1.5PN ~ The present paper will continue 
Papers I & II by investigating the gravitational radiation field of compact binaries, notably 
the total gravitational wave energy flux and orbital phase evolution, up to the same next- 
to-next-to- leading 3.5PN level. 

Including spin effects, and most importantly spin-orbit effects which are linear in spins, 
in the templates of gravitational waves emitted by compact binaries is of crucial importance 
for the accurate data analysis of the advanced ground-based as well as future space-based 
gravitational wave detectors. Astrophysical stellar-size black holes [3|-0] as well as super- 



massive black holes |8HlO| have spins, and the spins will affect the gravitational waves of 
black hole binaries through a modulation of their amplitude, phase and frequency (including 
the precession of the orbital plane in the case of non-aligned spins, see e.g. Refs. |lll-ll3|). 



The leading spin-orbit and spin-spin contributions in the equations of motion have been 

the next-to-leading corrections are also known both 



obtained using various methods |14|-|19 
for the spin-orbit 



20|-|24|] and spin-spin terms 



16, 
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: the 
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after a 



next-to- next-to-leading spin-orbit 
corrections have been derived in Refs. 28|, |29| and in Papers I & II. Concerning the radiation 
field of compact binaries, the leading spin-orbit and spin-spin terms are known 
next-to-leading spin-orbit terms at order 2.5PN ~ 1/c^ were first obtained in Ref 
previous attempt in |30l|: the 3PN ~ 1/c^ spin-orbit contribution including the tail integrals 
were computed in [3l|, after intermediate results at the same order (but including spin- 
spin terms) were given in 32| ; finally, the next-to-next-to- leading order contributions in the 
multipole moments and the energy flux, corresponding to 3.5PN ~ 1/c^, is the topic of the 
present paper and has never been addressed before. 

Following the previous investigations j2ll . 3l| we shall apply the so-called multipolar 
post-Newtonian approach to gravitational radiation. This approach has been extensively 



developped over the years (see [33|] for a review). It combines a mixed post-Minkowskian 



and multipolar expansion for the gravitational field in the exterior of a general source having 



compact support |34|-l37l|. with a matching to the post-Newtonian expansion of the inner 
field in the near zone of a post-Newtonian source [31 



42|. The gravitational waveform and 



various fluxes like the energy flux (or gravitational "luminosity" ) are expanded in a series of 
radiative multipole moments, that are then related to appropriate source-rooted multipole 
moments, expressed as integrals over the matter and gravitational fields in the source. The 
time derivatives of the multipole moments are performed using the equations of motion of 
the source, which must thus be known beforehand with the same accuracy as the one aimed 
for in the radiation field. 

Though the formalism can be applied to any post-Newtonian source, it does require a 
model for the source. In the case of compact binaries, the compact objects are described 
by point-like particles characterized only by their masses and their spins. The appropriate 
model in this context is an effective "pole-dipole" description, based on a stress-energy tensor 
made of a monopole or mass term involving delta functions, and a dipole or spin term made 
of gradients of delta functions. The pole-dipole description of spinning particles has been 



^ Hereafter these works will be referred to as Papers I & II respeetively. 
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developed by many authors j43|-|53j . The model is to be suplemented by an ultra-violet (UV) 
regularization in order to remove the infinite self-field of the point masses, see e.g. [s^, 55 |. 
In Papers I & II and the present paper, we use (and present justifications for that use at 
the aimed level of accuracy) the Hadamard "partie finie" regularization [56|, together with 
the Gel'Fand-Shilov prescription for distributional derivatives js^l, equivalent to Schwartz 
distributional derivatives |58|. Explicit expressions for the waveform and fiux can then be 
obtained in terms of the source's positions and velocities, and these can finally be converted 
into gauge-invariant quantities that are directly in use for building the gravitational wave 
templates. 

The plan of this paper is as follows. The Section [TTl is devoted to the general wave 
generation formalism and the Section HTTl to the application to spinning compact binaries. We 
give the required formulas in Sec. Ill Al for the various types of multipole moments (radiative, 
canonical and source), and in Sec. Ill Bl for the general post-Newtonian solution we employ. 
The complete results for the source multipole moments at next-to-next-to-leading spin-orbit 
level are presented in Sec. IIII A| while we obtain in Sec. IIII Bl the next-to-next-to-leading 
energy fiux and orbital phase evolution. We also present a numerical estimate of the new 
terms in Sec. IIII Bl In Appendix |A] we give alternative expressions of the source terms to be 
inserted into the source multipole moments. In Appendix [B] we check the agreement with 
the so-called boosted Kerr black hole limit, obtained when the mass and spin of one of the 
two black holes are set to be exactly zero. 



II. GRAVITATIONAL WAVE GENERATION FORMALISM 



A. Radiative and source multipole moments 



The gravitational waveform, generated by an isolated source described by a stress-energy 
tensor with compact support, is the transverse-tracefree (TT) projection of the metric de- 



viation, say hjj = {gij — 5ij) . It is defined in a suitable radiative coordinate system 
X'^ = (cT, X), at the leading-order 1/R when the distance R = |X| to the source tends 
to infinity, with the retarded time Tr = T — R/c being fixed. In radiative coordinates__the 
retarded time Tr asymptotically coincides with a null coordinate. The waveform reac 
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(2.1) 

The waveform is parametrized by two sets of symmetric and trace-free (STF) multipole 
moments, Ul of mass type and Vl of current type, which constitute the observables of the 



gravitational wave at infinity and are called the radiative moments [59| . They are functions 
of the retarded time Tr in the radiative coordinate system. Plugging Eq. (12. ip into the 



^ We denote by L = ii ■ ■ ■ 
from 1 to 3. Similarly L 



a multi-index composed of £ multipolar spatial indices ii, 



,ie rangmg 
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f_i and kL — 2 = kii ■ ■ ■ ig-2\ 



■ Ni^ is the product 



of £ spatial vectors Ni. In the case of summed- up (dummy) multi-indices i, we do not write the £ 
summations from 1 to 3 over their indices. The transverse-traceless (TT) projection operator is denoted 



v: 



TT 
ijkl 



'P^kVjl - \VijVki where Vi 



-NiNj is the projector orthogonal to the unit direction N = X/i? 



of the radiative coordinate system X'' = (cT, X). The quantity Sijk is the Levi-Civita antisymmetric 
symbol such that £123 = 1. The symmetric-trace-free (STF) projection is indicated using brackets or a 
hat. Thus Ul ~ Ul = t^(L) and Vl ~ Vl ~ V^l) for STF moments. We denote time derivatives with a 
superscript (n), and we indicate the symmetrizat|5)n operation with round parentheses. 



standard expression for the gravitational- wave energy flux we get 
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The radiative moments Ul and Vl are then related to some speciflc source-rooted multi- 
pole moments as follows. To implement the non-linearities in the propagation of the gravi- 
tational waves from the source to inflnity, we express them as some non-linear functionals, 
which can in principle be developed at any order, of some "canonical" moments Mj an d 5"^. 



When developped at the 1.5PN order they display the effect of tails and read [36|, |39 



Ul{Tr) 



2GM 



dtM 



(£+2) 



ln( 



t 



2GM 



Tr 



In 



V 2ro 
2ro 



o 



1 



(2.3a) 
(2.3b) 



where M is the mass monopole or Arnowitt-Deser-Misner (ADM) total mass. The quantities 
K,£ and 111 denote some numerical rational fractions and Tq is an arbitrary constant time scale; 
we shall not need any of these here. 

Next the canonical moments Ml and Sl themselves are given as some non-linear func- 
tionals of the "source" moments II and Jl, and also of four supplementary "gauge" moments 
Wl, Xl, Yl and Zl. In general, the canonical and source moments agree up to the 2.5PN 
order, namely 



Mr 



O 



1 



Sl = Jl + 



(2.4a) 
(2.4b) 



Since we address the computation of the spin-orbit 3.5PN contribution to the energy flux 
(12. 2p . we see that we only have to consider a possible spin-orbit term in the mass 
quadrupole (since spin contributions add at least a factor 1/c). The relation between Mij 
and lij is given by (see e.g. j60[): 
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(2.5) 



As was already noticed in Ref. [31|, the leading order spin contributions to W and lij both 
start at 0(c~^), so that {Mij)s = {Iij)s + 0{c~^). For our purposes, we can therefore ignore 
the distinction between canonical M^, Sl and source II, Jl moments. 

Finally the source multipole moments are deflned for a general post-Newtonian matter 
source for any multipolar order i > 2, and up to any post-Newtonian order. They are 
explicitly given by |4ipl 
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^ The brackets surrounding indices denote the STF projection; the STF product of £ spatial vectors is 
written as xl = X(^i^ ■ ■ ■ x,^^ = STF[a:;L]. 
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c^{i + + 2){2£ + 5) 

JL{t) = FP eab<^, j d=^X {r/r^f j dzi^5,XL-l>a S 



^m^L^g^l (x,t + zr/c) , (2.6a) 



2£ + 1 



■ 5^+1 XL_i>ac I (x, t + 2 r/c) , (2.6b) 



c2(£ + 2)(2£ + 3) 

The finite part operation FP in front represents an infra-red (IR) regularization defined 
by analytic continuation in a complex parameter 5, and involves tlie same arbitrary scale 
To = ctq as in Eqs. (12. 3p . wfiicti will be irrelevant for the present workQ 

The basic building "blocks" S, Sj and Sjj entering the latter formulas are evaluated at 
the position x and at time t + zr /c in a harmonic coordinate system (t, x) covering the 
source (where r = |x|). They are defined by 

S = - , Sj = , Sjj = T , (2-7) 

c 

together with r** = 6ijT^K Here r^'' denotes the post- Newtonian expansion of the total 
pseudo stress-energy tensor of the matter and gravitational fields, say 

r 1 

r^^'^PN W^"" + A^^'^ih) , (2.8) 

167rG J 

where T^*^ is the stress-energy tensor of the matter source, and A'^^{h) represents the gravita- 
tional source term which is given by a complicated non- linear, quadratic at least, functional 
of the field variable h^'^ and its first and second space-time derivatives. The pseudo-tensor 
appears in the right-hand-side of the Einstein field equations, when "relaxed" by the condi- 
tion of harmonic (or de Bonder) coordinatesJl The expressions (12.61) involve an intermediate 
integration over the variable z, with associated weighting function 

H^) ^ (1 - ^"Y ' 6e{z) = 1 . (2.9a) 

In practice the post-Newtonian expansion of the source moments (12. 6 p is performed by means 
of the formal infinite series 

+ 00 



E(x. /c) = g p^j^±|L_ 0V'>(x.<). (2.10) 



B. Explicit solution for the post-Newtonian metric 

To get explicit results at a given post-Newtonian order we need a solution of the relaxed 
Einstein field equations. As in Paper I, we parametrize an explicit solution by means of a 

^ This scale enters the relation between the retarded time in radiative eoordinates and the one in source- 
rooted harmonic coordinates: Trr ~ t ~ - — ^Si^ In ( — ) . 

^ The post-Newtonian expansion of the relaxed Einstein field equations takes the form Dh'^'^ = ^^^t^"^ , 
where □ = "qf^dp^ is the flat space-time d'Alembertian operator. Here /i'"^ = \/—g g^'^ — where g^'' 
is the inverse and g the determinant of the usual covariant metric g^^u] V^'^ is an auxiliary Minkowskian 
metric, 77'"^ = diag(— 1, 1, 1, 1). The harmonic-coordinate condition reads d^h^^ ~ 0. Note that the 
conservation of the pseudo tensor, d^jT^^ ~ 0, is the consequence of the harmonic-coordinate condition. 
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set of retarded potentials denoted by V, Vi, Wij, Ri, X, Zij. Here, contrarily to our previous 
work on the equations of motion, we will not need the highest-order potentials Yi and T that 
enter respectively goi at 0{c~'^) and goo at 0{c~^). All these potentials are "Newtonian" in 
the sense that they admit a finite non-zero limit when c ^ -|-oo. They enter the components 
of the usual covariant metric as follows: 



V' + - X + V,Vi + 



o 



1 



(2.11a) 
(2.11b) 



gij = Si 
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(2.11c) 



Equivalently, and more useful for the present work, they enter the "gothic" metric as: 
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(2.12a) 
(2.12b) 
(2.12c) 



Each of these potentials is a retarded solution of a flat space-time wave equation sourced by 
matter densities components and appropriate lower order potentials. The matter densities 
are defined from the components of the matter stress-energy tensor by 
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Denoting with n/j^S* the usual retarded fiat d'Alembertian integral, i. e. the retarded solution 
of 00 = ri'^'^dfj,u4> = -S", the latter potentials are defined by 



V = □;^^[-47rG(T] , 
= n^i[-47rGa,], 



AnG Van + WijdijV + 2VidtdiV + Vd^V 



+-{dtV)^-2diVjdjV, 
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Wi^ = [-AtxG (a,, - 5,,auK) - d,Vd,V] , 

- AnG V {a.j - bi^a^^) - ^df^iVdtVj) + diVkdjVk + dkVidkV, 



(2.14a) 
(2.14b) 



(2.14c) 

(2.14d) 
(2.14e) 
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-2d(iVkdkVj) - 5ijdkVm{dkVm - d^Vk) - -6,j{dtVf 



(2.14f) 



With the latter exphcit post-Newtonian solution in hands one obtains the basic buiding 
blocks (12.71) entering the source multipole moments (12. 6p as 
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(2.15a) 
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dtVdy - 2VjdtdjVi - Vd^Vi - 2djVdjRi + dtWijdjV - -VidjVdjV 



-2Vd,VAV - WjkdjkV. + d,W,AV, + dkW,,dkVj - dkWdkVi 
+2djVdiRj + 3VdtVdiV + VjdjVdiV + 2VdjVdiVj - duVjdiWjk 

+ 2' 
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(2.15b) 
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a., + — <: -U^A^ydkY + \d.yd,v 



"^{dtVY - dtVAV - \dkVAVk + \dkVidkV^ - dkVAV, 
dtVjdiV + dkV^d.Vk + dtVid.V + dkV^Vk - d,Vkd,Vk^ + O , (2.15c) 



where we have used the notation W = Wu and Z = Z^. Equivalent expressions for S, Sj 
and Sjj, which we have used for testing our calculations, are provided in Appendix Rl 
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III. APPLICATION TO SPINNING COMPACT BINARIES 



A. The pole-dipole effective formalism 



The pole-dipole formalism |43H53| is an effective description of point particles endowed 
with intrinsic (classical) angular momenta or spins, and moving in an arbitrary curved 
background — in practice the space-time generated by the particles themselves. The spins 
can take any orientation and magnitude, and in particular be close to extremal. In the 
present work we shall confine the formalism to terms linear in the spins. At that level 
the model can be used for describing black holes as well as ordinary compact bodies like 
neutrons stars. Indeed, the internal structure of the spinning bodies should appear only at 
the quadratic level in the spins, e.g. through the rotationally induced quadrupole moment. 

The stress-energy tensor of each of the particles is the sum of two terms, respectively 
built with a Dirac delta function and a gradient of a delta function, and integrated over the 
world line of the particle, according to: 

(3.1) 

Here the sum is over the two particles, r is the proper time measured along the world 
line of each particle given by the particle position ?/^(r); 5^^-' denotes the four- dimensional 
Dirac delta function; u'^ = dy^/ (cdr) is the four velocity of the particle normalized to 
UfiU'^ = —1; p'^ is its four linear momentum; and S^'^ denotes the antisymmetric tensor that 
represents the spin of the particle]^ Using a 3-1-1 space-time split, the particle's position and 
coordinate velocity are denoted y'^ = (ct,y(t)) and v'^{t) = (c, v(t)) (where = cu^/vP, 
with = [—gpaV^V^/c^Y^^"^), and we have 



dr V. 



^ ^/-gM c 



V-^(x,t) 



(3.2) 



where 5^^^ is the three-dimensional Dirac delta function and the spin tensor S^^{t) is con- 
sidered a function of coordinate time. 

Since the spin tensor S^'^ has six independent components, one must impose a supplemen- 
tary spin condition (SSC) in order to correctly describe the three independent components 



of the spin vector. Here we adopt the covariant SSC (Tulczyjew's condition |47|, |48| ) 

S'^V = . (3.3) 

It can be shown that, with the latter choice for the SSC, the mass defined by m?c^ = —p^p^ 
as well as the (four-dimensional) magnitude of the spin defined by S"^ = are 
conserved along the particle's trajectory. Furthermore the link between the four velocity 
and the four linear momentum is entirely specified. However, at linear order in the spins 
the linear momentum is simply proportional to the four velocity, 

pf" = mcu^' + 0{S^) , (3.4) 

^ In our convention the spin tensor and all spin variables have the dimension of an angular momentum times 
the speed of light c. This is useful for counting the effects of spins in the post-Newtonian approximation, 
as the spins appear to be formally "Newtonian" for maximally spinning particles. All the powers of 1/c 
are kept explicitly in our calculations, so that the next-to-next-to-leading spin-orbit terms we are looking 
for here will all carry in front a factor 
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so that the covariant SSC reduces to S^^u^ = 0{S^). Furthermore it can be shown that the 
spin precession equation reduces to the parallel transport equation 



DS 



dr 



(3.5) 



where D/dr is the proper time covariant derivative, while the four dimensional acceleration 
of the particle is given by 



2mc 



(3.6) 



This equation of motion is well known as the Mathisson-Papapetrou equation (43H46 



B. Multipole moments in the center- of- mass frame 

We work in all intermediate calculations with the spatial components and S2 of the 
two spin tensors (z,j = 1,2,3), eliminating the Oz components thanks to the covariant SSC 
fl3.3p . and then express all our results in terms of some spin vectors 5*^ and 6*2 . Like in Paper 
II we shall adopt a spin vector whose magnitude is constant in the usual three-dimensional 
Euclidean sense, namely 

5ijS\S{ = si , 5ijS\S2 = s\ , Si^2 = const . (3.7) 

We refer to Section IIA in Paper II for the precise construction, starting from the spatial 
components of the spin tensors, of spin vector variables 5*^ and 5*2 with constant magnitude, 
Eq. 03.71) . Such spin vectors therefore satisfy ordinary-looking precession equations, 

^ = e'^'' n{ ^ = e'^'' ni Si . (3.8) 

All the problem of the evolution of the spins SI and 5*2 in a binary system reduces to that 
of finding the ordinary precession vectors Q\ and fig. Those have been obtained for the two 
particles up to 3PN order in Section IIB of Paper II (for spin-orbit effects the precession 
vectors are independent of the spins). Now that we have defined the two constant magnitude 
spins SI and S^ it will be convenient to express the results in the center-of-mass frame by 
means of the variables [13] 

S' = Sl + S' = — S' S] . (3.9) 

7712 mi 

Here we denote the two individual masses by mi and m2, the total mass by m = mi + m2, 
and later we shall use the symmetric mass ratio u = mim2/m? (such that < z/ ^ 1/4), 
and the mass difference 5m = mi — m2. 

We compute the spin-orbit terms in the various source multipole moments which will 
enter the flux up to 3.5PN order in an arbitrary frame. At this stage, an interesting check 
described in Appendix [B] can be performed. Next, we adopt the frame of the center-of-mass 
(CM) defined in Paper II; notably the particle's trajectories y\ and 7/2 the CM frame are 
given in terms of the relative position = yl — 7/3 and relative velocity f* = dx^/dt = v\—vl 
by Eqs. (3.3) in Paper II. Finally the spin parts of the multipole moments are expressed 
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in terms of the conserved-magnitude spins and of the variables (13.90 . Our notation is e.g. 
(vS) = V ■ S for the ordinary Euchdean scalar product, e.g. (x x S)* = e'^^^x^T/' for the 
ordinary cross product, and e.g. {S,x,v) = S ■ (x x v) = e'^^^S'^x^v'^ for the mixed product. 
For each source multipole moment, we only present the linear in spin part, indicated by the 
underneath label S. We obtain, extending Ref. [2l| to next-to-next-to-leading order and 
writing all terms using our conserved spin variables. 
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We were able to perform a few technical tests on this calculation of the source multipole 
moments. We verified that the so-called "surface terms" can be computed either by a "bulk" 
integral over the entire three-dimensional space like other non-compact support terms, or by 
a surface integral extending on a sphere at spatial infinity. For such a test we have to use the 
alternative form of the blocks S, Sj and Sjj given in Eqs. (lAip : we refer to Section IVD of 
Ref. 6l| for a discussion of this type of terms. In addition we verified that certain quadratic 
non-compact support terms can be alternatively evaluated using some particular analytic 
kernels (denoted Yl, Sl and Tl in Ref. (s^). Another, more physical, test of the expressions 
of the multipole moments we obtain before going to the CM frame, is the agreement with 
the so-called "boosted Kerr black hole limit" as investigated in the Appendix [Bl 
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C. Flux and orbital phasing for circular orbits 



For the spin-orbit effects at the post-Newtonian level considered in the present paper we 
can neglect all the corrections 0{l/c^) in the relations between the canonical and source 
multipole moments, see Eqs. fl2.4p and fl2.5p . Furthermore the relations between the radiative 
and canonical moments, Eqs. fl2.3p . imply a spin-orbit contribution due to gravitational wave 
tails and arising at the 3PN order; we ignore this contribution here since it has already been 
computed in Ref. j3l|, and since the next-to-leading tail contribution would enter the result 
at 4PN order only. Finally, for our present purpose, we can replace all the radiative moments 
Ul and Vl by the corresponding source moments II and Jl up to the 3.5PN spin-orbit level. 
We can therefore use for the flux (12. 2p at that order the expression 
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The other terms do not contribute to the spin-orbit effect at the 3.5PN order. We insert 
the explicit results (I3.10p for the source multipole moments into Eq. (I3.1ip . we compute the 
time derivatives using systematically the equations of motion derived in Papers I & II, and 
we specialize the result to the case of quasi-circular orbits, again using the material from 
Papers I & II. 

It is useful to introduce an orthonormal moving triad {n. A, £} defined by n = x/r, 
£ = Lj<j/\L-^\ where = mi>:x. x v denotes the Newtonian orbital angular momentum, 
and X = £ X n. Then the spin-orbit contributions in the flux will depend only on the 
projections of the spins perpendicular to the orbital plane, namely Si = £-S and = € ■ S, 
where we recall that S and Xl are defined by Eqs. (13. 9p . Furthermore we denote the relevant 
post-Newtonian parameter for circular orbits by 



X 



(3.12) 



where u is the orbital frequency, related to the orbital separation r by Eq. (4.2) in Paper II. 
We are then left with the main result of the present work, namely the spin-orbit contribution 
to the flux up to order 3.5PN, as follows: 
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The tail-induced spin-orbit effect at 3PN order computed in Ref. [3l| has also been added, 
but we recall that we neglect spin-spin interactions. We have checked that this result is 
in complete agreement in the test-mass limit where — )■ with the result of black-hole 
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perturbation theory on a Kerr background obtained in Ref. 

To obtain the evolution of the orbital phase for quasi-circular orbits we shall apply the 
energy conservation balance equation relating the flux T to the energy E that is associated 
with the conservative part of the equations of motion: 



dE 
d^ 



(3.14) 



Note that the balance equation ( 13.14p is valid in average over a long radiation-reaction time 
scale ujjOj x~^/'^ = O(c^); thus short periodic variations at the orbital frequency u and at 
the spin precession frequencies Wprcc ^ xuj have been averaged out. In order to apply the 
balance equation fl3.14p we must ensure that the spins, or rather their projections Su and 
S2e (or equivalently Si and He), are secularly constant over the radiation-reaction time scale 
uj/u. 

This will be the case of the spin variables with conserved magnitude, as can be shown 
either explicitly at a given post-Newtonian order 6J| , or by the following structural argument 



valid at linear order in spins, extending the presentation of Ref. [31| . In the center-of-mass 
frame, the only vectors at our disposal, except for the spins, are n and v. Recalling that 
the spin vectors are pseudovectors regarding parity transformation, we see that the only 
way spin-orbit contributions can enter scalars such as the energy E or the flux J-" is through 
mixed products {n,v,Si) and (n, f, 6*2), i.e. through the components Su and S2e. Now, 
the same argument applies for the precession vectors r2i,2 introduced in Eqs. (13. 8p : they 
must be pseudovectors, and, at linear order in spin, they must only depend on n and v, 
so that we must have fli oc £ and fl2 oc £; this is explicitly seen for instance in Eq. (4.5) 
of Paper II. Now, the time derivative of the components along £ of the spins are given 
by dSu/dt = Si ■ [d£/dt + £ x fli] and idem for 2. The second term is zero, and since 
d£/dt = 0{S), we obtain that 5*1^ and 5*2^ are constants at linear order in the spins. This 
argument is valid at any post-Newtonian order and for general orbits, but is limited to 
spin-orbit terms. 

The conservative energy E has been obtained in Paper I and was reduced to circular 
orbits in Eq. (4.6) of Paper II. We recall here its expression: 
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Applying now the balance equation (I3.14p . in which we can assume by the previous argument 
that the spin projections 5*^ and are constant, one obtains the secular decrease of the 
orbital frequency as 
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Finally by a further integration we obtain the secular evolution of the orbital phase, or more 
precisely the so-called "carrier" phase defined by = J u dt, as 
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In the case of precessional binaries, for which the spins are not aligned or anti-aligned with 
the orbital angular momentum, the total phase $ is the sum of the latter carrier phase and 
the precessional correction arising from the precession of the orbital plane, $ = + 0prcc- 
The precessional correction 0pi.ec can be computed numerically [l3[ or analytically (see for 
instance Ref. [3l| for a computation at the IPN order). 

As a useful diagnosis to assess the importance of the latter spin effects, we have computed 
the number of accumulated gravitational-wave cycles between some minimal and maximal 



TABLE I. Spin-orbit contributions to the number of gravitational-wave cycles A/gw = (0max — 
0min)/'?'" accumulated from Wmin = vr x 10 Hz to Wmax = "^isco = (? I (6'^/^Gm) for binaries detectable 
by ground-based detectors LIGO and VIRGO. For each compact object we define the magnitude 
Xa and the orientation Ka of the spin by Sa = G Xa Sq and Ka = Sa • For comparison, we give 
all the non-spin contributions up to 3.5PN order; however we neglect all the spin-spin terms. 





1.4M0 + 1.4M0 


1OM0 + 1.4Mq 


lOJ\/0 + lOA/0 


Newtonian 
IPN 

1.5PN 
2PN 

2.5PN 
3PN 

3.5PN 


15952.6 
439.5 

-210.3 + 65.6kiXi + 65.6k2X2 
9.9 

-11.7 + 9.3K1X1 +9.3K2X2 
2.6 — 3.2kixi — 3.2K2X2 
-0.9 + 1.9K1X1 + 1.9K2X2 


3558.9 
212.4 

-180.9 + 114.0K1X1 + 11.7K2X2 
9.8 

-20.0 + 33.8K1X1 + 2.9K2X2 
2.3 — 13.2kixi — 1.3re2X2 
-1.8 + ll.lKixi + 0.8ft2X2 


598.8 
59.1 

-51.2 + le.OKixi + 16.0K2X2 
4.0 

-7.1 + 5.7kixi + 5.7K2X2 
2.2 - 2.6K1X1 - 2.6K2X2 
-0.8 + 1.7K1X1 + 1.7K2X2 
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frequencies, corresponding to the bandwidth of ground-based detectorsEI The results are 
given in Table [H They show that the 3.5PN spin-orbit terms computed in the present paper 
can be numerically larger, for spins close to maximal and for suitable orientations, than 
the non-spin 3PN or 3.5PN contributions. We thus conclude that they are still relevant 
to be included in the gravitational wave templates of LIGO/VIRGO/LISA detectors for an 
accurate extraction of the binary parameters. 
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Appendix A: Alternative expressions for S, Sj, Sj 



In this Appendix we provide alternative expressions for the building blocks S, S^, Hy, 
defined by Eq. (12.71) . which are useful for our practical computations. They are obtained 
by rewriting some products of derivatives so as to make Laplacians appear, for instance 
using 2diAdiB = A{AB) — AAB — BAA, and then using fl2.14p to replace the Laplacians 
of elementary potentials by source terms. In this process, new compact-supported terms 
appear, proportional to a, ai, aij. We find 
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^ Note however that the number of cycles of the carrier phase (j) docs not reflect the precession of the orbital 
plane, which has to be taken into account through $ = + 4'prcc- 
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+d,vdtW,, - 2V,dtd,v, + V,9,V - Ivd^v^ + ^ (^) ' ^^^^^ 



Gc^ \ 2 ' ' A ' 8 ^ " ' 2 



Such a rewriting allows an independent (and also faster and easier) calculation for all the 
"Laplacian" terms appearing in the multipole moments, for which the integrals regularized 
by means of the Finite Part operation take a simple form in terms of an angular average at 
infinity (see the discussion in Section IVD of Ref. (6l|). 



Appendix B: The boosted black hole limit 

In contrast to the case of the equations of motion for which we can perform several 
crucial verifications (see Paper I), there are not so many tests one can do in the case of the 
gravitational waveform and flux. In addition to the test-mass limit — of the flux f l3.13p 



which as we have seen perfectly recovers the result from Kerr black hole perturbations [63 



we can perform another physical test directly at the level of the multipole moments, before 



the reduction to the center-of-mass frame: the so-called boosted black hole (BBH) limit [66 

The BBH is obtained in the limiting case when we suppress one of the two black holes 
(say 2) by setting its mass m2 and spin S2 to be exactly zero into the general expressions 
of the multipole moments valid in an arbitrary frame, before going to the CM frame. What 
remains are then the multipole moments of a single Kerr black hole having mass mi and 
spin S'^-', and moving with constant velocity v\. In the BBH limit the multipole moments 
should agree with those of a single Kerr black hole moving with constant velocity, i.e. a 
black hole space-time on which a special Lorentz transformation or boost has been applied. 
The BBH test is interesting because it verifies (although only partially) the global Lorentz 
invariance of the multipole moments and the radiation field 



We start with the Kerr metric in harmonic coordinates. Since we are interested in spin- 
orbit effects we can work at linear order in the spin of the black hole. The "gothic" metric 
deviation H^^ of the black hole in the rest frame associated with some harmonic coordinate 
system = (cT, X) (thus satisfying d^H'^'^ = 0) reads 



= 1 - 'V 1:1/ . (Bla) 



GM 



H'' = ~-^ + 0(S-), (Bib) 



where M is the mass of the black hole and Si is the spin vector of the black hole in the rest 
frame. For convenience in this Appendix we define the spin vector directly from the spin 
tensor as Si = \eijkS^^. In the rest frame of the black hole we denote the spin tensor by S^^ 
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and the spin vector appearing in Eqs. (IBip is defined by Si = ^EijkS^^. We also denote the 
radial distance and unit direction in the rest frame hj R= |X| and = X^/R. 



Following Ref. 



we apply a global boost A'*j^(V) with constant velocity V = (V^). In 
our conventions we pose 



7: 



A^o(V)=A°,(V) = 7-, 
yiy. 



(B2a) 
(B2b) 

(B2c) 



with 7 = ( 1 



-1/2 



In the global frame defined by x'^ = h.^^X'^ the Kerr black hole 



metric will then be given by 



(B3) 



where (A ^)^^ = KJ^ denotes the inverse Lorentz transformation. The radial distance r and 
unit direction in the global frame = {ct, x) are related to the rest-frame counterparts 
R and A^^ by 



R = r 



r 



l + c^(f -1)(-) -2f(V^n)(-)+fi^ 



1/2 



n' - 7VM - 

r 



7^ V' 



(Vn) 



(B4a) 
(B4b) 



7 + 1 

See 6^; we denote the usual Euclidean scalar product by (Vn) = V^n^ . In addition the spin 
must also be transformed and we find that the spin vector Si in the rest frame is related to 
the spin vector Si in the global frame by 



Si — Si — 

7 + 1 c2 



(B5) 



Notice that the spin vector Si we are using here should rather be viewed as a covector, 
since it agrees with the spatial components of the covariant vector S"^ satisfying u'^Sfj, = 0, 
i.e. such that iSq = in the rest frame of the black hole. If one were to use instead the 
spin vector Si with conserved magnitude (as we did in all of this paper), one should have to 
apply a correction which is given at the 2PN order and for the BBH case by 



Si — Si 



(SV) 



8 c2 



Vi 



(B6) 



We now compute all the required multipole moments of the BBH by inserting the boosted 
Kerr metric (lBll) -( lB3l) into the general definitions of the source multipole moments. However 
this calculation is not straightforward starting from the defining expressions fl2.6p of the 



source multipole moments. Instead it was found in Ref. |66| that the best is to use some 
different expressions of the multipole moments, entirely given by surface integrals at spatial 
infinity; they have been derived in Eqs. (2.19) of 66|, see also (2.29) for the practical 



implementation. In the present paper we have used the same expressions of the multipole 
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moments and inserted there the BBH solution (lBl|) - (lB3p . Expanding the results at the 
required post-Newtonian order, and keeping only the spin parts of the multipole moments, 
we obtain: 



- (S X V)<,V,> 



4 
'3 



83 



Ji 



c 

t 



I ijk 

S 



+ ^y<^yJ>isv) 

= ^ (S X V)^,V,V, 



3 17V^ 

2 ^ 28 c2" + 336 



' 3^^ 
3 _ 5V 
1 4 c2 



6V^ 

7 c2 126 c4 
149 V^' 



J ijk — »S'<jVjVfc> 
S C 



k> 



o 



O C 



I ijkl 

s 



S 



^^(Sx V)<,V^,l^fcl^> + ofi 



5c3 



2 c 



(B7a) 



(B7b) 
(B7c) 
(B7d) 
(B7e) 
(B7f) 



These results are in perfect agreement with those obtained directly from our general com- 
putation of the multipole moments of black hole binaries, by setting m2 = and 5*21 = 0, 
and making the identifications mi = M, Su = Si, y\ = V^t and v\ = V^. 
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